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THE COMPUTATION OF THE PARTIAL VOLUMES 
OF THE COMPONENTS IN SOLUTION 

INTRODUCTION 


The partial volumes of the components in solution are important and are 
frequently used in tnermodynamic calculations. A number of methods have 
been developed for calculating these quantities [ 1] but the calculation of the 
partial volumes from the corresponding apparent volumes of the solution are 
convenient [2, 3] and appear to enjoy special favor because of a possible 
theoretical basis [ 4] . This is because of the empirical relation of Masson 
[ 5], who found that the apparent volume of the solute of electrolytes in aqueous 
solutions is a linear function of the square root of the volume concentration. 
Though this relation has found some theoretical confirmation for dilute solutions, 
it has often been found inadequate because the relation does not always extend 
to concentrated solutions. Modifications of Masson's relation are, therefore, 
frequently necessary, so that the essential problem of finding the partial vol- 
umes from a related quantity which can always be expressed analytically in the 
same form has not been achieved. It is not inappropriate, therefore, to describe 
the method given below for calculating the partial volumes since it is generally 
applicable to both aqueous and non-aqeuous solutions of both electrolytes and 
non-electrolytes. 


We have previously [ 6] derived the equation 


(P) (P) 


( p ) a (P) 

V s = X i° (V w ' 

(P) 


b) + bx t + x 2 ip 2 


(P) 


(i) 


where v , v and ip are, respectively, the specific volumes at pres- 

sure P, of the solution, water and the volume of the solute in solution. Also 

Xj and x 2 are the weight fractions of water and solute and a, the concentration 

of the solute, will be defined as either x 2 , the weight fraction of the solute, or 

by (x 2 /v ), the volume concentration (we will distinguish this by a'); b is 
s 

the constant in the equation 


P q - P = a log 


v (P >-b 

w 


( 2 ) 


, (P o> -b 


w 


which represents the variation of volumes of water with pressure. 



In this expression v^ P °^ and v^^ are the specific volumes of water at the 

w w 

pressure P and P, and the two parameters a and b may be determined by 

adapting the equation to the compressibility data of water. 1 Although equation 
(i) , and hence equation (2) , is expressed for water as solvent, the equations 
are generally applicable to both aqueous and non-aqueous solutions of electro- 
lytes or non-electrolytes. 

Equation (1) has been obtained by the combined use of equation (2) , 
which is a variation of Tait's equation [4, 6-9] and Tammann's hypothesis, 
to the effect that water in solution is compressed by an ionized solute. 

In particular, at atmospheric pressure, we will write for equation (1) 

v s = X 1 c “ ( v w ~ b ) + bx l + x 2 h (1') 

where for convenience we have omitted the superscript which designates pres- 
sure; v , v w and tp 2 will be used to indicate, respectively, the volumes of the 

solution, pure water and the volume of the solute at atmospheric pressure. We 
will derive our equations for the partial volumes for atmospheric pressure, 
the results holding approximately at higher pressures. 

Inequation (1’) only c and tp 2 are unknown (after a is appropriately 
defined and the value of b selected) and these may be found by writing this 
equation for each concentration (using density data for the solutions at atmo- 
spheric pressure) and determining by the method of least squares their best 
values. 2 


J The values of these parameters for water (and for a number of other 
solvents) have been determined with the initial pressure taken at atmo- 
spheric pressure; see Rosen [5] . b varies slowly as the pressure range is 
increased, but, as will be noted later, the value of b is not critical for our 
application. 

2 In reference [6] the parameters for a number of electrolytes in water 
were derived in this way from equation (1'). Equation (1) was then used 

( P) 

(after appropriate provisions were made for determining , the solute 

volumes at higher pressures) to compute the volumes of these salts in solution 
to 10 000 bars. 
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If we set 


ib . i = c® (v - b) + b 
ri w 

then equation (1’) becomes 

v s = X 1 + x 202 

where is the specific volume of water in solution (whose values depend on 
the concentration) and where ip 2 is the volume of the solute in solution and 
varies with pressure but is assumed to be constant for all concentrations. 

From the expression for in equation (3) it is evident that only a portion of 
the water, v^ - b, is compressed (since c< 1 for aqueous solutions of electro- 
lytes). Thus, b is the relatively incompressible portion of the volume of water 
and, because its value is not critically determined by equation (2) , neither 
is ip 2 , the solute volume in equation (1 T ). In addition, it is evident that the 
choice between the two definitions of a (or others that may be used) would 
determine the values of and ^ 2 . 3 However, the equations we will derive 
for the partial volumes, though different for either definition of the concen- 
tration used, give essentially the same numerical values for the partial 
volumes with different values for b and a. 


EQUATIONS FOR PARTIAL VOLUMES 


We rewrite equation (l r ) 

v = X! c“ (v - b) + b Xl + x 2 ^ 2 + 8, 

o W 

where the <5 is the difference between the observed and calculated volumes of 
the solution and is readily computed for each concentration once c and 0 2 
are found as previously indicated by the least squares method. 


3 For an interpretation of the volumes of the solute, 
from equation (i T ) for a particular choice of a andb, see 


ijj 2 , which result 
Rosen [ 9] . 
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The specific partial volumes of the components in a solution are given 
by the expressions 4 


▼l = v s ‘ x 2 (dv g /dx 2 ) 

v 2 = v + xj (dv /dx 2 ) 
s s 

where the subscripts 1 and 2 refer to the solvent and solute, respectively. 

Substituting equation (4) into these equations, we get the following 
relations for the partial volumes: 

^ - fa = -xjx 2 (v - b) (c“ln c) + A vj (6) 

W 

v 2 - fa = x 2 j (v - b) (c^ln c) + A v 2 (6') 

w 

and 

Vi - * -Xi (v - b) c a 'o!' (Vj/v ) Inc + Avj (7) 

w s 

v 2 - fa = (Xi 2 /x 2 ) (V - b) c“ a' (Vi/v s ) In c + Av 2 (7') 

where fa is given by equation (3) and the value of c corresponds to using for 
a and a ' the concentrations x 2 and (x 2 /v ). To these equations have been added 

the part resulting from the least square residuals, 5, and which we have 
designated by Avj and Av 2 ; so that by equations (5) and (5') can be written 


Avj = 

8 - x 2 (d<5/dx 2 ) 

(8) 

Av 2 = 

8 + xj (d6/dx 2 ) 

(8') 


The derivatives 5 inequations (8) and (8') are most readily obtained by 
numerical differentiation from a table of divided differences (of equal or unequal 
intervals of x 2 ) of the <5's as arguments of x 2 . In this table we are not limited 
to the arguments of x 2 used in fitting equation ( i') . This has several advan- 
tages in both convenience and accuracy. In the first place, the exactness in 
determining the parameters of equation (f) is not important, nor is it neces- 
sary to use all the density data available in obtaining this fit if the deviation 


( 5 ) 

(5') 


4 These define the partial volume per gram; the partial molal volume is 
obtained by multiplying these values by their molecular weight. Throughout 
this paper we use specific volumes instead of the corresponding molal volumes. 

5 <5 is of the order lxl 0 -4 (often, a fraction of this); however, dd/dx 2 is 
usually not negligible, especially at x 2 = 0. 
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curve is used. 6 Secondly, as new or revised data for the solution is obtained, 
only the contributions to the partial volumes given by equations (8) and (8') 
need be recomputed. 

When x 2 = 0, v< = v , but equations (6') and (7') (which becomes in- 
6 1 w 

determinate) quite readily give the equations 

(v 2 ) _ = ^ 2 + (v w - b) (In c) + (d6/dx 2 ) x2= 0 (9) 

x 2 “ 0 

and 

<v a ) = o = (V w - b) (l/v w ) In c+ (d6/dx 2 ) X2 = Q (9’) 

for a and respectively. At this concentration the partial volume is equal to 
the apparent volume and a stringent test for any interpolation equation is the 
limiting value obtained for this parameter. 

For x 2 = 1, equations (6 T ) and (7 T ) give v 2 , while equation (i f ) gives 

v = tp2* This limiting value for the partial volume has been identified as 
s 

the hypothetical maximum value of the apparent volume and associated with it 
is a "critical disruptive volume" corresponding to the upper limit to the con- 
centration of the solute in solution. For the alkali halides this comparison has 
been made with the limiting results derived from Masson’s empirical relation 
on apparent volumes [4, 9] . 

As we have previously indicated, the equations we have derived for the 
partial volumes also hold at higher pressures. In these equations all the 
volumes are expressed for the required pressure while c, a 9 a' and b are 
held constant. ^i(P) is given by equation (3), while fe(P) maybe obtained 

by a method previously described [4, 6] . 

In Tables I through III are shown the partial volumes of the components 
of aqueous solutions of ammonium nitrate. These, computed by the equations 
developed here, are compared with the carefully determined values given by 
Adams and Gibson [ 10] . For the method they use in determining these partial 
volumes see footnotes to Table I, Except at x 2 = 0, these values agree almost 
exactly at all concentrations. From the results in these tables it would also appear 


6 This is not peculiar to equation (1 ! ) alone; it is generally true when 
a deviation curve is used. However, this fact is helpful here as equation (1*) 
may be fitted by an approximate method more easily than by the exact least 
square method. 
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that our equations for the partial volumes of the components of the solution are 
not dependent on the value assumed for the parameter b, or whether the 
concentration of the solution is expressed by x 2 or x 2 /v g . 

At x 2 = 0, our values for the partial volume of ammonium nitrate 7 
differ from the value given by Adams and Gibson (0. 5932) , but out results 
are approximately what would be obtained if their partial volumes for the 
higher concentrations were extrapolated to x 2 = 0 ( 0.5981). 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, July 27, 1967 
129-04-01-00-62 


7 These are 0. 5976, 0. 5975 and 0. 5982 for the parameters shown in 
Tables I, H, HI, respectively. The slope of our deviation curve at x 2 = 0 
was computed from the data for 1, 5 and 10 percent solutions. 
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TABLE I. THE PARTIAL VOLUMES OF AMMONIUM NITRATE 
AT 25°C COMPUTED BY EQUATIO N (6') a 


*2 

b 

v s 

<5 C 

(v 2 ) C 

A” C 

Av 2 

v 2 

Eq. (6') 

- d 

v 2 

A & G 



10 4 A C 

0. 00 

1. 00293 

0. 00000 

0. 60337 

-0. 00573 

0. 59764 

0. 5932 

-44.4 

0. 15 

0. 94433 

-0. 00005 

0. 62011 

0. 00066 

0. 62077 

0. 6207 

-0.7 

0.30 

0. 88843 

0. 00001 

0. 63308 

0. 00015 

0. 63323 

0. 6332 

-0. 3 

0. 375 

0.86138 

0. 00003 

0. 63829 

0. 00002 

0. 63831 

0. 6383 

-0. 1 

0. 45 

0. 83489 

0.00001 

0. 64273 

-0. 00016 

0. 64257 

0. 6426 

+0. 3 

0. 525 

0. 80893 

0. 00000 

0. 64646 

-0. 00009 

0. 64637 

0. 6463 

-0.7 

0. 55 

0. 80040 

-0. 00003 

0. 64755 

-0. 00014 

0. 64741 

0. 6474 

-0. 1 


NOTES: 

a. The constants used in equations (1') and (6') are: 

b = 0. 8591, c = 0. 69215, ip 2 = °* 6563, a = x 2 . 

b. The data of Adams and Gibson [10]. 

c. The symbols used above are: 

6 is the difference between the observed 
specific volume of solution and the value 
calculated by equation (i 1 ) . 

(v 2 ) is used here to designate the partial 
volume calculated by equation (6 1 ) without the 
Av 2 term. 

Av 2 is given by equation ( 8') and is the correction 
to the partial volume caused by the deviations 6. 

A is the difference between the partial volumes 
found by Adams and Gibson and those computed 
by equation ( 6 T ) . 

d. The partial volumes computed by Adams and Gibson [10], 
They found the partial volumes by fitting two equations for the 
apparent volumes as functions of the weight fraction of salt and 
by using a deviation curve. 
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TABLE H. THE PARTIAL VOLUMES OF WATER IN 
AMMONIUM NITRATE SOLUTION AT 25°C 
COMPUTED BY EQUATION (7) & 


X 2 

b 

v s 

6° 

(vi)° 

A" c 

Avj 

Vl 

Eq. (7) 

-d 

Vl 

A&G 

10 4 A C 

0,00 

1. 00293 

0.00000 

1. 00293 

0. 00000 

1. 00293 

1. 00293 

0 

0.15 

0. 94433 

-0. 00020 

1. 00183 

-0. 00031 

1. 00152 

1. 00145 

-0.7 

0. 30 

0. 88843 

0.00001 

0. 99818 

-0. 00041 

0. 99777 

0. 99782 

0.5 

0. 375 

0. 86138 

0. 00008 

0. 99521 

0. 00002 

0. 99523 

0. 99523 

0 

0. 45 

0. 83489 

0.00003 

0. 99150 

0. 00071 

0. 99221 

0. 99222 

0.1 

0. 525 

0.80893 

-0.00015 

0. 98692 

0. 00171 

0.98863 

0. 98868 

0.5 


NOTES: 

a. The constants used in equations (i') and (7) are: 


b= 0.8591, c= 0.67362, il> 2 = 0. 66313, a* = (x 2 /v o ) 

s 

b. The data of Adams and Gibson [ 10] . 

c. Refer to the footnotes in Table I for the definition of the 
symbols (or their analogues) used here. 

d. See footnote d, Table I. 
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TABLE HI. THE PARTIAL VOLUMES OF WATER IN 
AMMONIUM NITRATE SOLUTION AT 25°C COMPUTED BY 
EQUATION (7) & 


*2 

b 

v s 

6° 

(vi)° 

a- c 

Avj 

Vi 

Eq. (7) 

_ d 

Vi 

A& G 

10 4 

0. 00 

1.00293 

0. 00000 

1. 00293 

0. 00000 

1. 00293 

1. 00293 


0. 15 

0. 94433 

-0. 00023 

1. 00188 

-0. 00036 

1.00152 

1. 00145 

-0.7 

0. 30 

0. 88843 

0. 00002 

0. 99827 

-0. 00043 

0. 99784 

0. 99782 

-0.2 

0. 375 

0. 86138 

0.00010 

0. 99532 

-0. 00014 

0. 99518 

0. 99523 

0.5 

0. 45 

0. 83489 

0. 00008 

0. 99146 

0. 00072 

0.99218 

0. 99222 

0.4 

0. 525 

0. 80893 

-0.00015 

0. 98668 

0. 00189 

0. 98857 

0.98868 

1.1 


NOTES: 


a. The constants used in equations (1') and (7) are: 

b = 0. 7750,' c = 0.77436, ip 2 = 0.6652, a'= (x^Vg) . 

b. The data of Adams and Gibson [ 10] . 

c. Refer to the footnotes in Table I for the definition of the 
symbols (or their analogues) used here. 

d. See footnote d, Table I. 
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